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Abstract 

We study a family of memory-based persistent random walks and we prove weak conver- 
gences after space-time rescaling. The limit processes are not only Brownian motions with 
drift. We have obtained a continuous but non-Markov process (Zt) which can be easely ex- 
pressed in terms of a counting process (Nt). In a particular case the counting process is a 
Poisson process, and (Zt) permits to represent the solution of the telegraph equation. We 
study in detail the Markov process {{Zt,Nt); t > 0). 

1 The setting of persistent random walks. 

1) The simplest way to present and define a persistent random walk with value in Z is to 
introduce the process of its increments (Yj, t G IN). In the classical symmetric random walk 
case, this process is just a sequence of independent random variables satisfying P(Yt = 1) = 
JP{Yt = —1) — I for any t > 0. Here we shall introduce some short range memory in these 
increments in order to create the persistence phenomenon. Namely (Yt) is a {—1, l}-valued 
Markov chain: the law of Yt+i given J^t = o-{Yo, Yi, . . . , Yt) depends only on the value of Yt. 
This dependence is represented by the transition probability n{x,y) — ]P{Yt+i = y\Yt = x) 
with ix,y) G {-1,1}': 

0<a<l, 0</3<l. 

The persistent random walk is the corresponding process of partial sums: 

t 

Xt = ^ Y, with Xo = Yo = 1 or - 1. (1.1) 

Let us discuss two particular cases: 

• li a + l3 = 1, then increments are independent and therefore the short range memory 
disappears. {Xt, t G IN) is a classical Bernoulli random walk. 

• The symmetric case a = {3 was historically suggested by Fiirth [7] and precisely defined 
by Taylor [T3]. Goldstein [S] developed the calculation of the random walk law and 
clarified the link between this process and the so-called telegraph equation. Some nice 
presentation of these results can be found in Weiss' book [17] and [l^. This partic- 
ular short memory process is often called either persistent or correlated random walk 
or Kac walks (see, for instance, [5]). An interesting presentation of different limiting 
distributions for this correlated random walk has been given by Renshaw and Henderson 

2) Recently, Vallois and Tapiero [15] studied the influence of the persistence phenomenon on 
the first and second moments of a counting process whose increments takes their values in 



j 1 — a a \ 
"^i P 1-/3 J 



{0, 1} instead of { — 1, 1}. They obtained some nearly linear behaviour for the expectation. 
Using the transformation y —> 2j/ — 1, it is easy to deduce that, in our setting, we have: 

IE_i[X,] :=IE[X,iXo = yo = -l] = ^(t + l)-7^^(l-p'+'). (1.2) 

TE+,[X,] := IE[^,lXo = Fo = +1] = ^ {t + 1) ^ (1 - (1.3) 

An application to insurance has been given in [16) . 

It is actually possible to determine the moment generating function (see Proposition 16.41 in 
Section 

$(A,t) =IE[A-^'], (AgK;). 

However it seems difficult to invert this transformation; i.e. to give the law of Xt- 
3) This leads us to investigate limit distributions. It is well-known that the correctly normal- 
ized symmetric random walk converges towards the Brownian motion. Let us define the time 
and space normalizations. Let ao and /3o denote two real numbers satisfying: 



< Qo < 1, 

Let Ax be a positive small parameter so that: 
< ao + coA^ < 1, 



< A) < 1. (1.4) 
0<A) + ciA, <1, (1.5) 



where cq and ci belong to R (see in subsection 16 . 21 the allowed range of parameters). 

Let (Yt, t £ IN) be a Markov chain whose transition probabilities are given by the matrix: 

A_/ 1 — ao — coAa; ao + coAx \ , . 

Let {Xt, t G IN) be the random walk associated with (Yt) (cf. (HH])). Define the normalized 
random walk {Z^, s G At IN) by the relation: 

zf=A,X,/A., (At>0, A, >0). (1.7) 

Set (Z^, s > 0) the continuous time process obtained by linear interpolation of (Z^). 
We introduce two essential parameters: 

po = 1 — ao — /3o (the asymmetry coefficient), (1.8) 

'70 =00- ao. (1.9) 

In this paper, we will aim at showing the existence of a normalization (i.e. to express At in 
terms of A^,) which depends on ao, /3o, so that (Z^) converges in distribution, as A^, — > 0. 
Our main results and the organization of the paper will be given in Section (2] 



2 The main results 
2.1 Case : po = 1 

Obviously po = 1 implies that ao = 0o ~ 0, and the transition probabilities matrix is given 
as 

A / 1 - coAx coAx \ , ^ 
" =[ cA. l-c,A. ) (-0,0 0). 

In order to describe the limiting process, we introduce a sequence of independent identically 
exponentially distributed random variables {e„,n > 1) with IE[e,i] — 1. We construct the 
following counting process: 

k>l 
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where 



1/co if k is odd 
1/ci otherwise. 



(2.2) 



Finally we define 



Z^o.ci ^ j (2.3) 







For simplicity of notations, in the symmetric case (i.e. Co = ci), Nf" (resp. Z^") will stand 
for Nt"'"" (resp. Zt"'""). The process (Z^") has been introduced by Stroock (in [13] p. 37). 
It is possible to show that if we rescale {Z^°), this process converges in distribution to the 
standard Brownian motion. This property has been widely generalized. For instance Bardina 
and Johs [1] have given weak approximation of the Brownian sheet from a Poisson process in 
the plane. 

Theorem 2.1. Let Ax = At and Yo — Xq = — 1. Then the interpolated persistent random 
walk (Zf , s > 0) converges in distribution, as Ax 0, to the process {—Z^"'''^, s > 0). 
In particular if Co ~ ci, then {N^°) is the Poisson process with parameter co- 
if Yo = Xo = 1 then the interpolated persistent random walk {Z^, s > 0) converges m 
distribution, as Ax — > 0, to the process {Z^^''^° , s > 0). 

Proof See Section H □ 

Next, in Section |3l we investigate the process (Z'j:''''^^ , N^"'"^^ ; t > 0). In particular we 
prove that it is Markov, we determine its semigroup and the law of {Z^"''^-'- , N^°''^^), t being 
fixed. This permits to prove, when co = ci, the well-known relation (cf. [18], [5], [8], [9]) 
between the solutions of the wave equation and the telegraph equation. For this reason the 
process {Z'^"''^^) will be called the integrated telegraph noise (ITN for short). 
We emphasize that our approach based on stochastic processes gives a better understanding 
of analytical properties. 

We will give in Section [5] below two extensions of Theorem 12. II to the cases where (Yt) is 

1) a Markov chain which takes its values in {yi, . . . , y^}, 

2) a Markov chain with order 2 and valued in { — 1, 1}. 

2.2 Case : po ^ 1 

In this case, the limit process is Markov. We shall prove two kind of convergence results. The 
first one corresponds to the law of large numbers and the second one looks like functional 
central limit theorem. 

Recall that {Z^, t > 0) is the linear interpolation of (Z^) and po (resp. 770) has been defined 
by ([L8} (resp. 

Theorem 2.2. 1) Suppose that rAt — Ax with r > 0. Then Z^ converges to the deterministic 

limit -lins- when Ax ^ 0. 
l-po 

2) Suppose that rAt = with r > 0, then the process (^^, t > 0) defined by 



converges in distribution to the process {^t, t > 0), as Ax 0, where 



- 2.(-^ + jr^)t+ (1 - 7^)^. (2.4) 

(Wt, t > 0) is a one-dimensional Brownian motion, t — (co + ci)/2 andT— (ci — co)/2. 

Proof See Section H □ 

Gruber and Schweizer have proved in [10] a weak convergence result for a large class of 
generalized correlated random walks. However these results and ours can be only compared 
in the case ao = /3o- 
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Note that ^ 

1 - „ = ^ ao = or /3o = 0. 

Suppose for instance that ao = 0. Then /3o,co > and 

Ct ^ Zt + ^= and = — — t. 

Obviously, the diffusion coefficient of (^J") can also cancel when po = —1. 
Since Po ~ —1 ao = /So = 1, then co, ci < and 

= and (,t = -rrt. 

This shows that, in the symmetric case (i.e. co = ci), we have = 0. This means that the 
normalization is not the right one since the limit is null. Changing the rescaling we can obtain 
a non-trivial limit. 

Proposition 2.3. Suppose ao = /3o = 1; co = ci < and rAt — AiJ with r > 0. 

The interpolated persistent walk {Z^ , t > 0) converges in law, as A^, 0, to {y/—rcoWt, t > 

0) where {Wt) is a standard Brownian motion. 

Proof. See subsection 16.31 □ 

2.3 Organization of the paper 

The third section presents few properties of the process [Z^"''^^ , i > 0) which has been defined 
by (|2.3|l . Theorem [211] will be proven in Section|31 Section [S] will be devoted to two extensions 
of Theorem 12. II In subsection 16 . II we determine the generating function of Xt (recall that Xt 
has been defined by (ll.lf) ). This is the main tool which permits to prove Theorem 12.21 and 
Proposition 12.31 (see subsections 16.21 and 16. 3p . 



3 Properties of the integrated telegraph noise 

The aim of this section is to study the two dimensional process {Z^°''^^, N^°''^''-; t > Q) in- 
troduced in (|2.2p and (|2.3|) . In the particular symmetric case co = ci, the study is simpler 
since the process {N^° , t > 0) is a Poisson process with rate co {JEj{N^°) = cot) and Nq° = 0. 
However we shall study the general case. 



First, we determine in Proposition 13. li the conditional density of Z'j:"''^^ given N^°''^^ — n. 
As a by product we obtain the distribution of Z^°''^'^ (see Proposition I3.3|l . Second, we prove 
in Proposition 13. SI that {Z^"''^^ , N^"''^'^ , t > 0) is Markov and we determine its semi-group. We 
conclude this section by showing that the solution of the telegraph equation can be expressed 
in terms of the associated wave equation and {Z'j:°''^°)t>o. For this reason, {Z^''''^''-)t>o will be 
called the integrated telegraph noise (ITN for short). Recall that: 

Co + Ci _ Ci - CO ,„ ^ , 

-=^— , (3-1) 

Proposition 3.1. 1) P(Ar=0''=i = 0) = e-*=o and given N^"'"^ = 0, we have Z^"'"^ = t. 

2) The counting process takes even values with probability: 

F(N:°'^' = 2k) = ^^Tk(k"-1)\ ""'^^ ^ / - + z)''e^^dz, (3.2) 

and the conditional distribution of Zl^''"^^ is given by 

P(Z,^«-^i Gdz|iV,=«'=i =2fc) = ^(t-z)^-nt + ^)'e^^l[_,.tj(z) (fc>l). (3.3) 

3) The counting process takes odd values with probability: 

p(^co.ci ^ 2fc + 1) = ''"a^fc+i^ll)? W = f + zfe~'dz, (3.4) 
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and the conditional distribution of Z'i:"''^^ is given by 

edz|iV,='"=i =2fc + l) = ^(t-^)^-(< + ^)'e^'l[-M]W ('=>0). (3.5) 

Corollary 3.2. In the particular symmetric case cq = ci, the conditional density function of 
given = n is the centered beta density, i.e. 

forn^2k, k€N' : f^t, z) = X2k ^ ' (3.6) 



with 



forn = 2k + l, keN: f„{t, z) = X2k+ ^^ ^ ^^j^, M-t,t]iz), (3.7) 

^ 1 (2fc + l)! . . 

X2fe+1 X2fc+2 22fc+iB(fc+ l,fc+ 1) 22'^+i(fc!)2 

is the beta function (first Euler function) : B{r,s) = ^^^^^q^^. 

Proof of Proposition \3.1\ Associated with ?i > and a bounded continuous function /, we 
define 

A„(/)=lE[/(Zr'"^^)l{^^=o.ci_}]. 

a) When ti = 0, we obtain 

Ao(/)=IE[/(Z,="-^l)i{,<;,,,,}]. 

If t < Aiei, then Z^"'"^ = t and 

Ao(/)=/(t)P(t<Aiei) = /(t)e-*^°. 

b) When n > 1, using (|2.1|l we obtain 

A„(/) = IE |^/(^t''°'''^)l{Aiei + ... + A„e„<t<Aiei + ... + A„ + ie„ + i}] ■ 

If AiBi + . . . + A„e,i < t < AiCi + . . . + A„+ie„+i then 



(-l)°dli+ / {-l)du+...+ / ( 

''■^1^1 -/ Aiei + ... + A„_ie„_i 



+ / (-l)"du. 

iAiei+... + A„e„ 

Hence 

Z^"'"' = Aiei - A2e2 + Xses + ... + (-l)""'A„e„ + (-l)"(t - Aiei - ... - A„e„). (3.8) 
c) Evaluation of A^kif), k > 1. 

We introduce two sequences of random variables associated with (e„): 

4fe = 62 + . . . + e2fe, Cfc = ei + . . . + e2fe_i, (fc > 1). (3.9) 
By l|3.8p . (|2.2p and (|3.9p we obtain the simpler expression 

A2fe(/) = IE \^f{t - '2^'k/ci)l{io/co+i%/ci<t<i'^/co+i%/ci + e2k + l/':o}] ■ 

Note that from our assumptions, ^fe, and e2k+i are independent r.v.'s, and are both 
gamma distributed with parameter k. Consequently: 

A2fe(/) = ((fc - 1)!)2 / exp{-co(<- y/co -a-/ci)}/(t-2a;/ci)e~"'~^a;*^"^y'°"^dxdi/ 



fc!(/c 



-jjj y /(f - 2x/ci)2;'°"^(t - a;/ci)'° exp - ijijc 



where Dt = R+ n {y/cQ + x/c-i < t}. Using the change of variable z — t — 2x/ci, we obtain 
X — ci , t — x/ci — and 

AMf) = — Jji-)[ [^) eM{ci~co)z/2}dz. (3.10) 

FinaUy pTTO)) and ((33J imply ((S^J and IfX^ . 

d) Evaluation of A2fe+i(/) for fc > 0. The arguments are similar to those presented in part c). 
On the event ^^+i/co+Cfe/ci < i < ^^+i/co +Cfc/ci +e2fe+2/ci, we have: Z^"'"' = 2Cfc+i/co-t; 
this implies 

A2fe+i(/) = E [l{^o^^/cQ+5=/ci<t}exp - ci(t- Cfc+i/c() -Cfc/ci))/(2Cfe+i/co - t)j . 

Since ^jl+i and ^£ are independent and gamma distributed with parameter fc + 1 (resp. fc), 
we get 



^2.+. (/) = |^e-(-+-)'/^ /(.) ( ^) (^) %xp { (c, - co)./2}d.. (3.11) 



This leads directly to (|3^ and (|33|) . □ 
Let us recall the definition of the modified Bessel functions: 

(^72)"+^'" 



m!r(!^ + m + 1) 

m>0 ^ ' 

Proposition 3.3. The distribution of Z^°''^^ is given by 

P(Z^=°''=i e dx) = e-'^'Stidx) + e-^' f{t,x)li^t,t]{x), (3.12) 

where 



f{t,x) ^ i[y '^°''^'^_^+ ^Ji(Vcoci(t^-x^)) +co/o(Vcoci(t2-x^))]e"^ (3.13) 

Remark 3.4. Let us focus our attention to the symmetric case Co = ci. We can introduce 
some randomization of the initial condition as follows: let e be a { — 1, l}-valued random 
variable, independent from the Poisson process , with p := P(e = 1) = 1 — P(e — —1). It 
is easy to deduce from (|3.12p that we have 

F{eZt"/t G dx) = (pSi{dx) + {1 ~ p)S-i{dx) + g{t,x)dxy-''°\ (3.14) 

with 

g{t,x) = ^{/o(cofyr^) + i±^^=il£A(cotN/r^)}l[_ia] W 
and 5i{dx) (resp. S-i{dx)) is the Dirac measure at 1 (resp. —1). 

In the particular case p — 1/2, x —> g{t,x) is an even function. G.H. Weiss filS^ p. 393) 
proved (|3.14p using an analytic method based on Fourier-Laplace transform. 

Proof of Proposition \3.3i The proof is a direct consequence of the expression of Proposition 
13.11 Indeed, for each bounded continuous function tp we denote 



A = nAZt°'"')] = + J2 + E ^2fc+i(^) = ^(t)e-=°' + Ae + A„, 

fe>l fc>0 

where A„(^) = IE[v3(^r'''')l{jv,'^o.ci^„j]. Using (|321) and ((33} we get 

'^'y ip{z)Se{z)e'^''dz, 



Ae = e" 
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with 

(coCi)" (t- Z\k-T^ (t + Z\k 



2 ^ A:!(fc- 1)! 



1 , [tTz^ 1 /\/coCi(t2 - 22)N2fc + l 



' fc>0 ^ ' 



^ k\(k + 



2 

For the odd indexes, by (|3.4p and p.Sp we get 

Ao = e~^*y ifi{z)So{z)e^''dz, 

with 



fe>0 



□ 



Proposition 3.5. (Z^ "'''^ iV^'''''''^ ; t>0) is aRx '!N-valued Markov process. 

2) Lets>0 and n > 0. Conditionally on ZJo.<=i ^ ^ ^yf'-'^i = n, (^{Z^l'^^ , N^°-^''^), t > 0^ 

is distributed as 

+ /,'( — 1)^" du,n + Nf°''^-'-^ , t > when n is even, 
(^(^x- J*{-l)'^^^'''°du,n + Nt^-''°^, t>0j otherwise. 



Remark 3.6. Note that Propositions \3.5\ and \3. 1\ permit to determine the semigroup of 
^{Z^o,ci^^c„,c,-^^^ > o) I.e. FiZ^"'"^ e dec, iVt"«'"i = nlZ^"'"! = y, N^"'"' = m) where 
t > s, n > m and y £ [— s, s] . 

Proof of Provosition \3.5\ Let t > s>Q. Using (|2.3p we get 



0-1 p 
JO 



where iV^ = Ar^^-;;! - N",^'"^ ,u>0. 

Note that (iV^; u > 0) = {N^"'"^; u>Q) if iV^^o.^i g 2IN and (iV^; u> 0) = (iVS''''°; « > 0) 



if Af^O'^^i £ 21N+1. This shows Proposition 1331 □ 
Next, we determine (in Proposition 13.81 below) the Laplace transform of the r.v. Z^°''^''- . It 



is possible to use the distribution of Z'^"''^^ (cf Proposition 13. 3p . but this method has the 
disadvantage of leading to heavy calculations. We develop here a method which uses the fact 
that {Zl°''^'^; s > 0) is a stochastic process given by (|2.3p . The key tool is Lemma [3.71 below. 
Roughly speaking Lemma 13.71 gives the generator of the Markov process {Z^°''^^ ,N^''''^^). 
Lemma 13.71 is an important ingredient in the proof of Proposition 13. 1 11 besides. 

Lemma 3.7. Let _F : R x IN — > R denote a hounded and continuous function such that 
z — > F{z,n) is of class for all n. Then 

^^nF{Z:°''\N:°''')] =Ie[ ^{Z^O'^\N:'>'^'){-lf^'"^'] 

+ IE [ (^F(Z,"°^"i,iV,"°'"i +1) -F(Z,"«'"i,iVt"0'"i)) 

^ ('^i^{iVt"-"ie2]N+i} + colfAT^o.^igaiN})]- (3.15) 

Proof. Let us denote by A{t) = ]E[F{Z'j:"''^^ , N^°''^^)]. In order to compute the t-derivative 
we shall decompose the increment oi t A(^) in a- sum of two terms: 

A{t + h)-A{t) 

r ~ Jjh + Oh, 
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with 

Since F{-,n) is continuously difFerentiable with respect to the variable z and t —> Zl°''^^ is 
diflerentiable (cf (|2.3|l 'l. using the change of variable formula we obtain 



Therefore 

lim Bh^TE [^(^-o,<=i jv,'=«-=i)(-l)<'"'''l. (3.16) 
In order to study the limit of Ch, we consider two cases: Nt"'"^ £ 2 IN and N^"'"^ e 21N+1: 

+ i IE [(f(ZI°-^^ , Nt-'^- + Nl°-^- ) - FiZ^O'^^ , N:°-^^ )) Ij^co.e, J , 
where Nh = iVt+h' - N^"'"' . 



According to Proposition 13.51 conditionally on Z^"''^^ and N^"''^'^ £ 2 IN (resp. N^"''^^ G 
21N+1), Nh is distributed as AT^o-'^i (j.ggp iV^i-'^o). Note that Proposition \3A] implies that 
p(jY^o,ci > 2) = o{h) and 

Consequently 

hm = ciE [(F(Z,^«-\iV,^«-i + 1) - F(Zr"\^r'^^))l{<o-ie2^+i}] 

+ CO IE , iV,^«-i + 1) - F{Z:«'^^ , N^-'^^)) ^j] . (3.17) 

Then, (^1^ and (PTT)! clearly imply Lemma □ 
Let us introduce the two quantities: 

Le(t) =IE [e-'*^'""'l^^co.eig2^J andio(t) =IE [e-'*^'""'l{^co.cig2^^^j], (t > 0, ^ G R). 

(3.18) 

Since \Zf°''^^ \ < t, then Le(t) and Lo{t) are well defined for any ^ G R. Note that — > -Le(t) 
(resp. Lo{t)) is a Laplace transform. We have mentioned the t-dependency only because 
it will play an important role in our proof of Proposition 13.81 below. 

Proposition 3.8. Let L^it) and Lo{t) be defined by (|3.18p . Then 

L,{t) = (^(-^ + -) sinh(tv^) + VI cosh{tV£)'j e""', (3.19) 

Lo{t) = smh{tV£)e-^\ (3.20) 
]E[e-MZ?-"i ] = -L ^(_^ + ^) sinh(tv^) + VS cosh(fyf )j e""*, (3.21) 
where £ = /i^ — 2r/i + . 

Proof. Applying Lemma [3.71 with the particular function F{z,n) = e~'^^l{„g2iN}, we have: 
^ie(t) = -MlE[e--"'"(-l)""'"l{^co,.,^^^j] 

+ F^e ' (^l{iv"0'"is2]N+i} - l{]Vj"0'"ie2iN}J 
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We deduce 
Similarly 

d 



^L^t) = -(^i + C(,)Le(t) + ciLo(t). 



+ -^['^ * (^l{JVj''0'''l e2]N} ~ l{JVj''0' = l g2]N +1} j 



We get therefore 
To sum up 



^Lo W = (m - ci)Lo(t) + coLe(t). 



\ holt) J V '=0 J \ Loit) 

We deduce the expressions of Leit) and Lo{t): 

=a+e^+' + a_.e^-' = 6+e^+* + 6_e^-*, (3.22) 

where A± = — r ± ^ fi^ — 2r^ + = — r ± ^/S. 

The constants a± and b± are evaluated with the initial conditions: 

Le(0) = F(iVo"°'"i e 2IN) = 1, Lo(0) = P(Aro"°'"i G 21N+1) = 0, 

^(0) = -(^t + Co)ie(O) + CiL„(0) = - Co, 
or 

^(0) = - ci)L„(0) + coL.(O) = CO. 
at 

We obtain 

a+ = — r=(— /i + r + Vf) and a_ = — -=(/i — r + Vf), (3.23) 
2v5 2vf 

fe+ = and fe_ = %= (3.24) 

2^ 2\/£ 

Using dX^l}, (p:^ and (PIM)) . Proposition ED follows. □ 



It is easy to deduce two direct consequences of Proposition 13.81 First, taking /i = we 
obtain W{Nl°^''^ G 2W) and W{N^°'''^ G 2IN+1). Second, taking the expectation in we 
get the mean of Z'^"''^^ . 

Corollary 3.9. We have: 

p^^co,ci g 2IN) = i|^-sinh(rt) +rcosh(rt)]e"^*, 
p(jV^co,ci £ 21N+1) = — sinh(Tt)e-"*, 



and 



lE[Zr''^^] = -t+^(l-e- 



Remark 3.10. The Laplace transform with respect to the time variable can also he explicitly 
computed. We define F{ii,s) = Jjj°° e"** IE[e^^^*'' ^]dt. Integrating p.2ip with respect to dt 
we get 



2^/FV ^s-VS + r 2V£^ ^s + VS + T 

{VE + t){s + V£ + t) + {V£ + fi- t){s -VI + t) 

2V£{{s + r)-^ - £) 
2sV£+^TV£-2fi^/I _ s + 2t-h 
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In the symmetric case, £ equals jj,'^ -\~ Cq, then 

s + 2co — 
+ 2sco — j-i' 



(3-25) 



Let (Zt) he the symmetrization of {Z^°) which is defined by an initial randomization: 



Zt = eZt\ t>Q, 



where e is independent of Z^° and P(e — ±1) — 1/2. 
Relation p.25[) implies 



f 

Jo 



-^^^]dt- ^ + ^"° 



/ s2 + 2sco — fl'^ 

This identity has been obtained by Weiss in il8^ . 

Let us now present a link between the ITN process and the telegraph equation in the 
particular symmetric case cq — ci = c > 0. Recall that (Nf) is a Poisson process with 
parameter c. 

Let / : R — > R be a function of class whose first and second derivatives are bounded. We 
define 

u{x,t) ^ -^f{x + at) + f{x - at)^ x € K, t>0. 
Then (cf [5]) m is the unique solution of the wave equation 

'dW d^' 

u{x,0)^f{x), -^{x,0) = 0. 
Proposition 3.11. The function 

w{x,t) =M^u(x, J (-l)^'ds^j, (xeR,t>0) 

is the solution of the telegraph equation (TE) 

w{x,0) = f{x), ^(s,0)=0. 

This result can be proved using asymptotic analysis applied to difference equation asso- 
ciated with the persistent random walk [8] or using Fourier transforms [18]. Here we shall 
present a new proof. 

Proof of Provosition Iff, jil Applying twice Lemma 13.71 to {z,n) u{x,z) and {z,n) 
^{x,z){-l)" we obtain; 

^ix,t)^^[^{x, I i-ir=ds](~l) 

and 

=|2 



^ix,t) = IE [^{-,1 i-lf^ ds\] ^ 2c TR 
Since u solves the wave equation we have 



The function w is actually the solution of the telegraph equation. It is easy to prove that w 
satisfies the boundary conditions. □ 



Let us note that Proposition 13.111 can be extended to the asymmetric case co 7^ ci . In 
this general case the telegraph equation is replaced by a linear system of partial differential 
equations. 
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Remark 3.12. 1) In fS^, an extension of Provosition [37Tl\ has been proved. Let A be the 
generator of a strongly continuous group of bounded linear operators on a Banach space. If w 
is the unique solution of this abstract "wave equation": 

^=A'w- w{;0) = f, ^{.,Q)=Ag {f,geV{A)) 

then u{x,t) — IE j^w ^3;, (— 1)^= ds^ j solves the abstract "telegraph equation": 

_=^,_2c-, n(,0)=/, ^(•,0) = A5. 

8) In the same vein as f^, Enriquez fS^ has introduced processes with jumps to represent 
solutions of some linear differential equations and biharmonic equations in the presence of a 
potential term. Moreover useful references are given in 
3) It is easy to deduce from Lemma\3. 7| that the functions 



pi 

,(x,t)=]E\u(x, J (-l)'^=°'''ds)l{jvco,cig2N}]. ixeR,t>0) 



Wo{x,t) ='!E^u(^x, (-1)^"° 'ds^l^^oo,cig2]N+i}], {x&R,t>0) 
are solutions of the general telegraph system (TS) 



= (CoCl — CojTOe + (CQCl — Cj )uio — iCQ + d 



^ " ^ " • ' " ' " " " dt ' dx'^ ' 

= (CoCi - Co) We + (C()Ci - Ci )TOo - 2ci + a -^-^ , 

w^{x,Qi) = f{x), Wa{x,Q)=Q ^^[x,Q) = ~cof{x) ^^{x,0) = cof{x) 



4 Convergence of the persistent walk to the ITN 

Suppose po = 1. The aim of this section is to prove the convergence of the interpolated 
persistent random walk towards the generalized integrated telegraph noise (ITN) i.e. Theorem 
12.11 Let us start with preliminary results. 

First, let us recall that {Xn,n £ IN) is the persistent random walk starting in defined by 
the increments process {¥„, n £ IN) (see Section [T|) with transition probabilities 

A _ / 1 - coAt coAt 
^\ ciAt 1-ciAt 

Let (Tfc; fc > 1) be the sign changes sequence of times : 

J Ti = inf{t > 1 -.Yt^Yo} 

\ n+i = mi{t >n -.Yj^Yt,}; k>l. 



(4.1) 



We put To = and 

Ak=Tk-Tk-i k>l (4.2) 
Let Nt be the number of times over [0, t] so that the sign of {Y„) changes: 

Nt = J2kT,<t) (4.3) 

The definition of Nt implies that: 

Nt = k<i^Tk<t< Tfc+i 
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We suppose in this subsection that Yq = —1. 
We deduce from the identities above: 

k 

= + -Tfc + 1) where fc = iVt. (4.4) 

By (|4.2p we obtain: 

Tk = Ai + ... + Ak k>l. (4.5) 

Hence tlie equations (|4.3|) . (|4.4p and (|4.5|) permit to emphasize the bijective correspondence 
between (X„; n G IN) and (A^; fc G IN). 

We introduce the normalization of {X„; n G IN) given by p.7p with A^; = At: 

(s/AtGiN). (4.6) 
Let us define: 

N^ = J2^i^UA.A,<s} ->0. (4.7) 
J>1 

Let us note that 

Nf" = iV,/A, if s/At G IN . 
That permits to extend the definition of to any s > by setting 

k 

2f = ^(-l)''(At^) + (-l)'=+^(s-Atrfc + At) k = N^. (4.8) 

j=i 

Obviously = if s/At G IN. 

In order to study the asymptotic behaviour of (Z^) as At 0, we shall first prove the 
convergence in distribution of {AtAj)j>i and (Nf^)s>o- 

We recall that some random variable ^ is exponentially distributed with parameter A > if 
its density is given by j e~^^'^l^^:^oy. 

Lemma 4.1. The random variables (Ak) are independent and AtA2k (resp. AtA2k+i) con- 
verges in distribution, as At 0, to the exponential law with parameter ^ (resp. 

Proof. Since (Yn) is a Markov chain, then the {Ak) are independent. First let us study the 
convergence in distribution of the sequence AtA2k. We use the Laplace transform of AtA2k'. 
(fiifM) = IE[e~'"^''*2'=], /i > 0. Since A2k is geometrically distributed with parameter ciAt, we 
obtain 

oo 

V.(m) = ^e-^^"(l-CiAO^'-'ciAt 
i=i 

ciAt ciAt ci 



e"^' - (1 - ciAt) {^i + ci)At + o{At) 11 + ci 



-o{At) (4.9) 



The function ip{j-i) converges for any /i > to the Laplace transform of some exponential law 
with parameter c^^. This proves the convergence in distribution of AtA2k. Concerning A2k-i 
the arguments are similar. □ 

Let us recall that the counting process {N^"''^^ , t > 0) has been defined through the se- 
quence of jumps (e„;n > 1) via (|2.ip . and (e„; ?i > 1) are i.i.d. and exponentially distributed. 

Lemma 4.2. Let s > 0, fc > 1 and $fe : R'^ — > R 6e a bounded continuous function. Then 

1) ^lim^F(iVi^ = 0) = P(iV,"°'"i = 0) 

2) \\m.W.[^k{AtAi,AtA2, . . . , At4fe)l{jvA^fe}] = IE[$fc(Aiei, A2e2, . . . ,\kek)l[i^-o.<=i^k}\^ "^^ere 



Xk has been defined by (|2.2p . 

Proof. 1) Statement 1) follows from: 



p(ivf = 0) = P(iVu/Ad = 0) = P(ri > [s/Atj ) 

= F{Ai> [s/At\ ) = P(AMi > At Ls/AtJ ) 
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where [aj denotes the integer part of a. 

2) Set k>l. The event {Nf" = k} can be decomposed as foUows: 

r fe > r k+i 

{N^ = fc} = I A, ^ ^ < s| n I A, ^ A, > 

This identity imply existence of a bounded Borel function ■ R''^^ R so that 
$fc(AMi,...,Atylfe)l{ivA=fe} 

= $fe(AMl, . . . , AtAfc)l{At Ej=i A,<.}l{At Ej+i A,>.} 

= •(/'fc(AtAi, AM2, . . . , At^fe+i). 
Since $fc is continuous, the discontinuity points of tpk are included in: 

k fe + 1 

1^: G R'^+' : = sj U |a; G : I] 2;, = s|. 



U : 



By Lemma l4.ll {AtAi, . . . , AtAk+i) converges in distribution towards (Aiei, . . . , Xk+iBk+i) 
as At — » 0. Since the Lebesgue measure of U is null, the limit law does not charge U. We can 
conclude evoking for instance Theorem 14 p. 247 in [3]). □ 

Let us formulate a straightforward generalization of Lemma 14.21 

Lemma 4.3. Let n G IN, (fci, . . . , k„) G IN" such that ki < k2 < ■ . ■ < k„ and {si, . . . , s„) G 
R" mth Si < S2 < ■ ■ ■ < Sn. Let $ : R*" ^ R 6e a bounded and continuous function. Then 

^lim^lE[$(Aiyli,...,AtylfcJl{jvA^,^_..._^A 

= ]E[<l'(Aiei,...,Afc„efe„)l{jv|o.-i=fc^,...,jv|o,ci^j^^j] (4.10) 

Proposition 4.4. The random variable converges m distribution towards —Z^"''^^, for 
any s > 0, as At ^ 0. 



Proof. Let / : R — > R be a continuous function which is bounded by M. Identities (|4.8p and 
gSl) imply that IE[/(Jf )] = ^^=0 ^^(fc), with 

k k 

E^{k) = IE [/( ^(-1)^ AM, + (-1)'^+' - At E A, + At)) 1{^A=,}] 
j=i 3=1 



Applying Lemma 14.21 and p.8|l . we obtain for any fc > 0, 

k k 

^lini^£A(fc) = ij[/(^(-irA,e, + (-1)'+' (« - ^ ^^^0) ^{^^"■'''-^>] 

Moreover since / is bounded by AI, we get 

\EA{k)\ < MP(iVi^ = fc). 

Suppose that fc > 1. Then, using the Markov inequality and the independence property of 
the random sequence {A„, n > 0), we obtain 

k k+1 

P(iVf = fc) = p(At^A, <s< At^ylj 

j=i j=i 

< P(^At^Aj <s) =P(^exp|-At^^| > 



< e^IE[exp-AtX^^] ^e^Hv'.^l) 



i=i j=i 
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where (^j(^) = ]E[e-'"^"*^]. Since is a Markov chain starting at = — 1, for any j > 1, 
^2j-i (resp. A2j) is geometrically distributed with parameter coAt (resp. ciAt). According 
to (|4.9|l we get 

'^^^■(1) = eA._i + ciA. ^ A, + c,A, = < 
By the same way, we have: 

¥'2,-l(l)<^<l. 

1 + Co 

As a result, there exists < r < 1 so that 

F{N^ = fc) < e''r^ (4.11) 
We are now allowed to apply the dominated convergence theorem: 

hm IE[/(^f )] = J2 i5A(fc) = ElE[/(-^.''"'^)l{ivfo-i=,}] = IE[/(-^r«'^^)]. 



□ 

Proposition 4.5. For anj/ (si, . . . , s„) G R" siic/i i/iai si < S2 < ■ • ■ < Sn, fie random vector 
{Z^^, . . . , Z^^) converges in distribution to (— Zf^''^^ , . . . , — ), as At tends to 0. 

Proof. We follow the approach developed in the proof of Proposition 14.41 Let / : ^ R be 
a bounded and continuous function. We have: 



E 



[f{Z^^,...,Zt)]= J2 EA{ki,...,k„), 



where the sum is extended to (fci, . . . , G E^" so that fci < fc2 < . . . < fc„ and 

_Ba(/ci, . . . = IE [/(^fj, . . . , Z^)l{jvA^fej_..._jvA^=fe^}j. 
Identity (|4.8|l implies the existence of a bounded continuous function tpn '■ R'°" ^ R so that 

EAiki, . . . , fc„) = IE [i/;„(AMi, . . . , Atylfc„)l{ArA =fe,,...,ivA . 
Applying Lemma 14.31 we get 

lim -EA(fci, • . . , fc„) = lE['i/'n(Aiei, . . . , Afc„efe„)lrjv'^o.<:i^j, pf^o.^i^f. J. 

According to the definition of the process Z^"'^^ , we may deduce : 

^iim^i5A(fci, . . . = nfi-Zsr\- ■ . ,~^j°'^^)i{^jo.^i=fe^,....^jo.=i=,„j]. 

In order to obtain that 

hm n.f{zt ,---,zt)]- n.fi~z^,r' > • • • , -zz-^' )], 

Af -^0 

it suffices (cf the proof of Proposition 14. 4|) to prove that 



sup |£;A(fci, . . . , fc„)| < oo. 

Since / is bounded, 

|£A(fci,...,fc„)| < MF{Nfl =fc„) 
Using moreover (|4.11|l we get 

J2 |£A(fci,...,fc,0| < Afe""^(fc„)"-V'=" <oo 

ki,...,kji kn 

since r < 1. □ 
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sup \Z^-Zf\<5. 



We are now able to complete the proof of Theorem 12. II Since {Z^) and {Z'^°''^'^) are both 
contmuous processes, the convergence of the process {Z^) to the process (— Z^"''^^) will be 
proved as soon as the following measure tension criterium (cf Theorem 8.3 p. 56 in [2]) holds 
: for all e > and 770, there exists some constants 5 £]0, 1[ and 71 > such that 

i P ( sup \Z^ - Zf\ >e)< 770, for any At < ^i. (4.12) 

Since (Z^ , s > 0) is the interpolated persistent random walk, its slope is always equal to 1 
or —1. Hence we obtain for any (m, s) G R^, 

\Z^-Zf\ < \u-s\. 

Consequently 

s <u<s + (5 

By choosing 5 = e/2 we get the tension criterium and so the convergence of the process (Z^) 
to the process {-Z^O'"^). 

5 Two extensions of Theorem 12.11 

First of all, the extensions presented in this section concerns the regime A^ = At. 

5.1 The case when (Yt) takes k values. 

Let us introduce our parameters. Let fc > 2, yi, j/^ denote fc real numbers, and (c(i, j); 1 < 
hj < fc) a matrix so that 

c(i, j) > for any i / j, c(i,i)=0, ^c(i,/)>OVi. (5.1) 

We directly consider the asymptotic regime. Let (Yt) be a {yi, . . . , j/fc}-valued Markov chain, 
with transition probability matrix: 

= \ / . ^ . . (5.2) 

1- E,=ic(j,0 At 7 = J, 



where At > is supposed to be small so that 

k 



c(i,j)At<l, (^c(i,0)At<l 



Similarly to the case k = 2 and yi = — 1, 172 = 1, we are interested in the linear interpolation 
(Zf ; s > 0) of the process (Zf ; s > 0) defined by (tL7)) . 

Theorem 5.1. Suppose Yo = yi- Then {Z^; s > 0) converges m distribution, as At 0, 
to the process ^ Rgds; t > where (Ra) is a {yi, . . . ,yk}-valued continuous-time Markov 
chain starting at level yi, whose dynamic is the following: {Rt) stays on level yi an ex- 
ponential time with parameter 1/^ X^iLi c(i, 0^ "'^d, jumps to yji (j' 7^ j) with probability 

c(j,i')/(Etic(j,o). 

Remark 5.2. In the case fc = 2, yi = -1 and y2 = 1, then ((-1)^*"°'"'; t > 0) (cf (pTTIl j 

may be chosen as a realization of {Rt ) when it starts at Ro — —1. 

Proof of Theorem 15. Jl We proceed as in the proof of Theorem 12.11 developed in Section |4l 
Let (T„)„>i be the sequence of stopping times defined by (|4.1|l . Then; 

^ _ y^{t+l) 0<t<Tl 
y»Ti + Yt, (t-Ti + l) Ti < t < T2. 
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Recall that {Z^; s/At G IN) has been defined by (|4.5|l . From the relations above, it is easy 
to deduce: 

y,(s + At) 0<s<AtTi 

y,(AtTi) + Yt, (s - AtTi + At) AtTi < s < AtT2. 
Let us determine the limit distribution of (AtTi, VVi) as At — > 0. Set 

y^(A,i)=IE[e-^^'^il{y^^=,^}], \>0, j^i. 
Proceeding as in the proof of Lemma l4.ll we obtain: 



1-1- 



Using standard analysis, we deduce that {AtT\,YTi) converges in distribution as Af ^ to 
(e'l, L^i) where: 



IE |e ^'=il{y^=,}| = 



c(j, j) 



A + Y. 'i=i I) 

As a result, e'l and Ui are independent, e'l is exponentially distributed with parameter 
l/E!LiC(i,0 and 

P(C/,.,-) = _£(!ii^. 

Ei=ic(i,0 

Using the approach developed in Section [l] we can prove Theorem 15.11 The details are left 
to the reader. □ 

5.2 The case when (1^) is a Markov chain of order 2. 

Let (Yt) be a Markov chain with order 2. For simplicity we suppose that it takes its values in 
{ — 1, 1}. Obviously {Yt,Yt+i)t>o is a Markov chain with state space 

i5 = {(-1,-1), (-1,1), (1,-1), (1,1)}. 

Let TT^ be the transition probability matrix: 



A 



-coAt 


Co At 














1 — Po 


Po 


Pi 


1 — Pi 














ciAt 


1 - ciAt 



(5.3) 



where At, Co, ci, po,pi > and coAt, ciAt,po,pi < 1. 
Let us introduce: 



= Z 7Z ^^^-TTZ T, c[ = avi, i = 0,1. (5.4) 



l-(l-po)(l-pi) 

Recall that (Z^) and (Zf") have been defined by (|4J)) . resp. (g^, {Nt"'"^) is the counting 
process defined by (|2.ip . and 

^="•=1 = fi-lf'"'"^ du, t > 0. 
Jo 

Theorem 5.3. 1) Suppose that Yo = Vi = -1 (resp. Yo ^ Yi = 1) then (Z^; s > 0) 

converges m distribution, as At -> 0, to {-Zs"' ^ s > 0) (resp. [Zs^' °; s>0)). 

2) Suppose Yo = 1 and Yi = —1 (resp. Yo — —1, Yi = 1) then {Z^; s > 0) converges in 

distribution, as At 0, to 

((6 - 1) ^ Vl)^"°"'' du + ej\-lf^^^''''du- s > O) 

where e is independent from [Nu"'^^), {Nu^''^") and 

P(e = 0) = 1 - P(e =1) =vi {resp. P(e = 1) = 1 - P(e = 0) = vq). 
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Remark 5.4. 1) Note that {Yt)teT!si ts a Markov chain if and only if 1 — cqAj — pi and 
1 — CiAt =po- If we replace formally po (resp. pi) by 1 — ciAt (resp. 1 — coAt) in (|5.4|) and 
take the limit At 0, we obtain Vi = pi and c'i — Ci. We recover Theorem \2.1\ 
2) The fact that {Yt) is a Markov chain with order 2 does not modify drastically the limit. 
The limit process can be expressed in terms of processes of the type {Zf''^; s > 0). 

Proof of Theorem \5.y\ 1) We only consider the case Yq = Yi = 1. Let us define Ti, T2 and 
T3 as follows: 

Ti = inf{t >l,Yt = -1}, T2 = mi{t > Ti + 1, Yt = Yt-i}, T3 = mi{t > T2 + I, Yt / Yt^}. 
Using the definition (cf (HHI) of (Xt) we easely obtain: 



Xt = 



t + 1 < t < Ti 

Ti + Xt Ti<t<T2 



where Xt equals either —1 or 0. 
Moreover, when T2 < t < T3, we have: 



Xt 



Ti - 2 - (t - T2) if T2 - Ti is odd 
Ti + 1 + (t - T2) otherwise. 



Zf 



According to p.7p . we can deduce: 

' s + At < s < AtTi 

AtTi + AtX,/A, AtTi < s < AtT2 

AtTi - 2 At - {s - AtTa) AtTa < s < AtTs, Yt^ = -1 

AtTi + At + s- AtTa AfTa < s < AtTg, Yt^ = 1 

(note that Ta — Ti is odd if and only if Yt2 = — !)• 

2) a) Proceeding as in the proof of Theorem 12.11 we can prove that AtTi converges in 
distribution, as At 0, to e'l, where e'l is exponentially distributed with parameter 1/ci. 

Then (Zf ; < s < AjTi) (s; s< ei), as A* 0. 

b) The distribution of Ta — Ti does not depend on At. Moreover |X. | < 1, then the limit of 
the length of the interval [AtTi, AtTa] is null. We have 



p(yT, = -1) = ^((i-Pi)(i-Po))Vi = ^^1- 



c) Using the strong Markov property, we easely show that (^i^AtTj ; < s < At(T3 — Ti)) — > 
(e'l+yTiS; < s < Ca), as At — > 0, where {e'i,YTi) (resp. (e'i,ea)) are independent r.v.'s and 
conditionally on Yt2 = 1 (resp. Yt2 ~ —1) Ca is exponentially distributed with parameter 
1/ci (resp. 1/co). 

d) Let us summarize the former analysis. We have proved that {Z^; s > 0) (/^ Rudu, s > 
0), where (-Ru) is a continuous-time Markov chain which takes its values in { — 1, 1} and Ro = 1. 
Moreover the dynamic of (-Ru) is the following: (-Ru) stays in 1 (resp. —1) an exponential 
time with parameter 1/ci (resp. I/cq) and moves to —1 (resp. 1) with probability vi (resp. 
Vo)- Note that (-Ru) is allowed to stay in the same site. It is classical (cf [12]) to prove that 

{Ru)u>o == (^u^''^°)u>o where Cq and c'l are defined by (|5.4|l . □ 

6 Convergence of the persistent random walk to- 
wards the Brownian motion with drift 



In subsection l6.1l below we determine the generating function of Xt, where Xt is the persistent 
random walk defined by p.l|l . This allows to prove Theorem 12.21 and Proposition | 
subsections 16.21 16.31 
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6.1 The moment generating function of Xt 

Let us recall that the increments process (Yt, t G N) is a Markov chain valued in the state 
space iJ = { — 1, 1}. Its transition probability is given by 

/l — a a \ „ 
7r= 0<a<l, 0</3<l. 



f3 1-/3 

The persistent random walk {Xt, t G N) is defined by the partial sum: 

t 

Xt = ^ Fi with Xo = Ko = 1 or - 1. 

Lemma 6.1. Let us define the functions at and bt: 

atij)=-P{Xt=j,Yt = -l) and bt{j) = SPiXt = j,Yt = 1). (6.1) 

Then, 

at+i{j) = (1 - a)atij + 1) + (3bt{j + 1) (6.2) 
bt+i{j) = aat{j - 1) + (1 - f3)btU - !)■ (6-3) 
Proof. Using the Markov property of (Yi) we have: 

at+iU) = F(Xt+i=j,Yt+^ = -l,Yt = ^l) + nXt+i=j,Yt+r = -l,Yt = l) 

= P(Xt + Yt+i = -1, Ft = -1) + TP{Xt = j + 1, Ft+i = -1, Yt = 1) 
= (l-Q)at(j + l)+/3&t(j + l). 

The second recursive formula involving {bt{j)) can be obtained similarly. □ 

Let us define the moment generating function $(A, t) — 1E[A^'], (A > 0). We decompose 
$(A,f) as 

$(A,t) = $-(A,t) + $+(A,t), (6.4) 

with 

<I.-(A,t) = IE[A^n{y,=_i}], $+(A,t)=IE[A^n{y,=i}]. (6.5) 

Lemma 6.2. 1) ^-(A.O) = ^ P(Fo = -1) and <E'+(A,0) = AP(Fo = 1)- 

2) The moment generating function verifies the following induction equations: 

$_(A,i + l) = i^ $_(A,t) + ^ $+(A,t) (6.6) 
$+(A,t+l) =aA<l>_(A,t) + (l-/3)A<l>+(A,t) (6.7) 

Proof. Definition (|6.ip implies that 

cI>_(A,f) = ^A^at(j)= E ^'"'(^■)- 

jez j=-i-i 

Hence, 

<I._(A,t + l) = EA^at+i(i) = (l-a)EA^at(i + l)+/3EA^&t(j + l) 

j 3 3 

= (1 - a)i E A^+iat(j + 1) + f E + 1) 

J 3 

= ^^$-(A,t) + ^$+(A,t). 

The proof of (|6.7|l is similar. □ 
Lemma 6.3. Let f{\,t) be equal to either $_(A,t) or $+(A,t), then 

/(A, t + 2) - + (1 - /3)a) /(A, t + 1) + (1 - a - /3)/(A, f) = 0. (6.8) 
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Proof. By (|6.6|l . we get 



$+(A,t) = {$_(A,t + 1) - i^<J>_(A,t)}-. (6.9) 



Replacing f by f + 1 in (|6.9|) . we obtain 

$+(A,t+l) = |ci,_(A,t + 2)- i^$_(A,f + l)}|. (6.10) 

Using successively (|6.7p . (|6.10|l and (|6.9p . we have: 

aA$_(A,t) = $+(A,t + l)- (l-/3)A$+(A,t) 

= ^|cE,_(A,t + 2)-i^<l>_(A,t + l) 

^"^A^|$_(A,f + l)-i-^ci,_(A,t) 



Finally 



<I._(A,t 



+ 2) - (^i^ + (1 - /3)A^ $_(A,t + 1) + ((1 - a)(l - /3) - a/3)$„(A,t) = 0. 



The proof concerning /(A, t) = "l>_(A,t) is similar and is left to the reader. □ 

In order to obtain the explicit form of $-(A,t) and $+(A,f) in terms of A and t, it suffices 
to compute the roots i9_ and i9+ of the following polynomial 

^iy^ + (l-/3)A^ ji»+l-Q-/3 = (6.11) 

Its discriminant equals 

(iy^ + (l-/3)A)'-4(l-a-/3). (6.12) 

It is clear that 

^ = (i^ + (l-/3)A + 2^)(i^ + (l-/3)A-2Vp) 



A 

1 /I - a 
A 



(-j^ + (1 - P)\ + 2Vp) ((1 - /?)>' - 2VpA + 1 - a) . (6.13) 



Since the discriminant of A ^ (1 — l3)X^ — 2^fp\ + 1 — a is equal to — 4a/3 then V > for any 

A > 0. 

Consequently the roots of (|6.1H) are: 

^?± = i(^^ + (l-/?)A±^/p). (6.14) 
We deduce the following result. 

Proposition 6.4. 1) The moment generating function <I>(A,t) satisfies 

$(A,t) = a+i9+ + (6.15) 

with 

1-q + A(Aq-»9-.) 1 -fv V ^ 

a+ — r= ana a~ — — — a+ if Ao = lo = — 1 

A^vO A 

and 

(1 - /3)A" + /3 - Ai?- . . fv V ^ 

a+ — ;= and a- — \ — a+ if Xo = Yo = 1- 

Vv 
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Proof. Suppose that Xq = Yq = —1. Let us first determine the values of the generating 
function at time t — and t = 1: 

$(A,0) = '1>+(A,0) + $-(A,0) = jF{Yq = -l) + AP(yo = 1) = i = a+ + a_ 



Moreover, using (|6.6p and (|6.7|l with t — 0, we get 

$(A,1) = $+(A,l) + $_(A,l) = (^i^ +qA)$_(A,0) = +a = a+i9+ +a_t?- 



It is clear that Lemma lOl and <I>(A, t) = ^^.(A, + $-(A, t) implies that <1>(A, t) satisfies (|6.8p . 
Then (|6.15p follows by standard arguments. The second case Xq — Yq — 1 can be proved in 
a similar way. □ 

6.2 Proof of Theorem [HI 

We keep the notations given in Section [T] Let qq and /3o be two real numbers in [0, 1]. Let 
Aj; be a small space parameter so that: 

0<ao + coA^<l, < f3o + ciA^ < I, 

where cq and ci belong to R. 

Note that ao > (resp. /3o > 0) implies that ao + cqAx > (resp. /3o + ciAa; > 0) when A^; 
is small enough. If ao < 1 (resp. /3o < 1), similarly ao + cqA^ < 1 (resp. /3o + ciAx < 1) 
as soon as Ax is small. In the case ao = 1 (resp. /3o = 1) co (resp. ci) has to be chosen in 
]-c»,0]. 

We assume that the coefficients of the transition probability matrix tv'^ of the Markov chain 
(Yt) satisfy: 

a = ao+coA^, /3 = /3o + ciA^ (6.16) 

i.e. TT^ is given by (|1.6|l . {Xt) is defined by (|I.1|I and (Z^) is the normalized persistent 
random walk: 

= A,X,/A, , (A* > 0, A, > 0, s G At IN). 
{Zi^; t > 0) denotes the linear interpolation of (Z^). 

Recall that po = 1 — ao — /3o and 770 = /3o — ao- Note that po ^ 1 ao + /3o 7^ 
Proposition 6.5. Let po 7^ I, 

1) if rAt = Ax with r > then Z^ converges towards the deterministic limit — 
tends to 0. 

2) if rAt = A^ with r > 0, Z^ + converges in distribution to the Gaussian law 
with mean _ 



and variance 

2 r(l + Pq) ( r]l 



a — 



po (I-po)^- 



1 - po V (1 - /5a)2/ 
where 

c = Co + ci and c = ci — cq. (6.19) 

Proof. We shall prove the statement under the condition Xq = Yo = — 1- If ^0 ~ Yo = +1, 
the limit is obtained by changing the sign and replacing co (resp. ci) by ci (resp. Co). 
1) Let <I>(A,t) be the generating function associated with Xt. In order to determine the limit 
distribution of Z^, let us introduce: 

^{li,t)^E.i[e-^^^], (6.20) 
where ]E_i denotes the expectation when Yg = —1. Observe that 

0(p,t) =$(e-^^^-^) =E_i[e-^^^^'('/^')], (6.21) 

when t/At e IN. 

According to Proposition 16.41 when t/At £ IN, (j){p,t) can be expressed in terms of a+, a_ 
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and \pV. 

First let us study the asymptotic expansion of the discriminant D as A^; 0. It is convenient 
to set: 

I = coAj; and I = ciA^;. (6.22) 
Applying (|6.12p with a = ao + 5 and /? = /3o + 5 we have: 

D = ((1 - Qo - 5)e^^^ + (1 - /3o - 5)e"^^-) ' - 4(1 - qo - A) - ^ - <5). 
By (fOSj) we get 

O =( (2-ao-/3o)+A,(^^(A)-ao)-co-ci) (6.23) 



(^^ (2 - Qo - /3o) + m(ci - co)) + o(A^))' 
- 4(^l~ao-/?o- A^co + ci)) (6.24) 
It is clear that D admits the following asymptotic expansion, as A^, 0: 

7? = Ao + AiA, + A^^l + o(A^) 
It is usefuU to note that ao and /3o can be expressed in terms of 770 and po: 

ao = ^ and /Jo = ^ . 

Let us compute Aq, yli and A2 using standard analysis: 

^0 = (2 - QO - /3o)' - 4(1 - QO - /3o) = Qo + /3o + 2ao/3o = (qo + /3o)' = (1 - po)'. 

.4i = 2(2 - Qo - /3o) (^/i(/3o - Qo) - (co + ci)) + 4(co + ci) 

= 2^(2 - Qo - A))(/3o - Qo) - 4(co + ci) + 2(qo + /3o)(co + ci) + 4(co + ci) 

= 2|^t(2 - QO - A))(/3o - Qo) + (qo + A))(co + ci)| 

= 2(^7,0(1 + po) + c(l-po)). (6.25) 



2(2-ao -/3o)(^y(2-ao -/3o) +^(ci - co)) + (^p(/3o - Qo) - (co + ci))^ 
= ^^((2-ao-/3o)' + (/3o-ao)^) 

+2/1 (^(2 - QO - A))(ci - Co) - (A) - ao)(co + ci)^ + (c,, + cxf 
= 2At^((Qo-l)' + (/3o-l)') +4/t(^(l- A))ci - (l-Qo)co) +(co + ci)2 
= At^(ry,^ + (l + po)') + 2Ai(^(l + po)c-r?oc) (6.26) 
Under the condition po 7^ 1, we have 



Hence 
with 



= (1 - Po)Jl + (^A. + ^^Ai + o(Ai) 



= So + Bi A, + B2A^ + o(A^) 
Bo = 1 - po, 



Bi = 



1 ^1 



2 1 — po 1 — po 

7?o(l+Po) , 

1 - Po 



— i — |/i(2 - QO - /3o)(/?o - qo) + (qo + /3o)(co + ci)| 
1 - Po I- J 
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_ 1 1 Al 

As a result, B2 is a second order polynomial function with respect to the /i- variable: 

B2 = B22 + M^21 + S20- 

Identities (fOS)) . (|626]) and (fOT)) imply: 

(2c(l-po))'_ 
^ 2(1 -po) " 8(1 - por ^ ° 

^2((l + po)c-r,oc) l 877oc(l-po)(l+po) -l + Po 2»?oc 
i>2i = : — = c- 



1-/50 8 {l~P0f 1-/90 (l-/90)2 

lri + {l+p,f 14%2(l + po)^ 1 ('7o' + (l + Po)^)(l-po)^-%^l + Po)^ 

i>22 = X" ~ 



2 1-po 8 (l-po)3 2 (1-po)^ 

(1 + Po)^ _ 2ripo 
2(1 -po) (l-po)3- 

Consequently 

?7o(l + Po) , \. , / 2/(l + Pof ^VoPo 



PR 1 , /^ '7o(i + Poj , \ . , r 2 /(J^ - 



Po ' 7 ^ ' r V2(l-po) (l-po)3 

+p(^^-7T^)}a^+o(A^). (6.28) 
V 1 — po (1-po)'' /J 

2) The first order development suffices to determine the limit of (/!>(p, f) as A^: 0. Indeed 

A. 

po 

From (|6.14p and (|6.16p we can easely deduce 



Then 



= I-P0 + — ^|wo(1 + Pd)+c(1-po)|+o(A,). (6.29) 
1 - po J 

id (|6.16[l we can easely deduce 
^(1 + po) + ^(m. - c) ± l{l - PO + A.( ^-'°/^+;°) + c) } + o(A.). 



,9+ = l + A,-^+o(A,) and ^_ = po - A + c) + o(A,). (6.30) 
1 — po Vl — po / 

Let t' = Li^J At. Since = ^t'^ + (f - tOAa^yu/Ad+i and jy^] < 1, then 

\<i>{p,t)-4>{p,^)\<Ci\r.t.t, (6.31) 

where C is a constant which only depends on p. 
Recall that identity (|6.20p and Proposition 16.41 lead to 

<^(p, t') = a+^^'J^' + a-{i''J^' (6.32) 

where 

(1 - a)e^''^^ + a - j^-e**^^ 

a_|. = ^ -= and a- = e'' ^ — a+. (6.33) 



It is obvious that (|6.33|) and (|6.30|) imply: limA^^o a+ = 1 and limA^j^o a- = 0. 
Since limAt^o ^- = Po and —1 < po < 1 then 

lim a-'&^'l^* = 0. (6.34) 

Ai^At^O 

Consequently, the second term in (|6.32|) tends to 0. It is important to note that the initial 
condition Xq — Yq — —1 disappears. Let us study the first term in the right hand side of 
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(|6.32|l . Note that limAa:-.o ~ 1, then if A^^ is small enough, we can take the logarithm of 
From (|6.30|l a straightforward calculation gives 

log^+ = A,-^+o(A,) 
1 - Po 

Choosing rAt — and using (|6.32p . (|6.34|l and (|6.31l) . we obtain the following limit: 

hm (ji{fi,t) =exp{- }. 

Ax^O 1 — Po 

Since the convergence holds for any £ R, we can conclude (cf Theorem 3 in [4]) that 

lim ]E_i[exp(mZ(^)l = exp | — ^^''^o^ \ ^ f^j. g^^y u £^1. 
Ax^O L 1 — On J 



Thus Zt converges in distribution, as A^, 0, to the Dirac measure at — 
3) Next, we consider the convergence of the process 

^A ^ ^A , tr|o^/r 



(l-po)x/A7 
Hence we define 

i>{ti,t) =IE_i[e~''^* ] =e (i-''o)v^,^(^,t). 

To determine the limit of tp{fi,t) as At, A^, 0, from (|6.32|l and (|6.34p we may deduce that 
it suffices to compute the second order development of the root Using (|6.14|) and (|6.28p 
we get: 

Al ffi^l + po) 



-(1 + Po) + —{^J■vo-c) + —[ ^ +Mcj 



_j_ 1 - Po ^ A^ / ^?7o(f +po) ^ ^ 



+ 2 



2 2 V i-po 



As a result 



1?+ = 1 + A^- — + A^(— I —--r, m 71 vt) )+o(A^). (6.35) 

1 — Po V 2 VI — po (1 — po)'^/ VI — po (l-po)-'// 

We take rAt = A^. Then 

J;».*(M) = („,<,r^xp{-i^r^J-}) 

f Af^yoi 1 rt ^ „ 1 
= hm exp < - h -pr log } . 

It is straightforward to deduce 

2 2 

^lim '(/'(P) — 6xp I — nifi + ~^~| (6.36) 



with 

-c r;oC 



1 - po (1 - po) 



(6.37) 



.^ = !ll±M(l--^)t. (6.38) 

1 — po V (1 — poyj 

4) Since (|6.36p holds for any /i G R, this implies that converges in distribution, as A^ ^ 0, 
to the Gaussian distribution with mean m and variance a^. (see Theorem 3 in [4]) □ 
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Proposition 6.6. Assume that po ^ 1 and rAt = (A^;)'^. Let us denote the process 
defined by 



(l-po)x/A7 

Then (^^ , , . . . , ) converges in distribution, as Aj, — > 0, towards {^i^ ^J"^ ) where 

is given by 



.1-po (I-Po) 
(Wt, t > is the one- dimensional Brownian motion starting at 0. 

Proof. The proof is only presented in tlie case n — 2. For simplicity let s — ti < t2 = t. 
We are interested in the limit of the random vector Let us then compute the two 

dimensional Fourier transform 

*'^(m,A) =IE-i [e'^f^^'-^^'e'^^^], (A,mGK). 

Since the process {Xt,Yt) is Markovian, we obtain 

*^(m,A) = ^-l[e''^^-^^]■^E-l[l^Y(s/A,)^^l}e^^^] 

+ E+i [e'''^^"-] IE_i [l{y(VA,)=+i}e'^^^] , 

when s/At and t/At belongs to IN. 

Note that - < Aa;At when m' = At. Consequently 

+ IE+1 [e'^«^^-='] IE_i [l^Y(s'/A,)=+l}e^^^'], (s' = [s/AtJA*, t' = lt/At\At). 
According to Proposition 12.21 

lim E-i\e"'^^-^'] = lim £+1 [e'^«^^-='l ^ e^^"""-^ "'^^'-'^ 

A^^O I J A^-»0 L J 

where m and are defined by (|6.37p . resp. (|6.38|) . Then we can deduce: 
lim *^(/i,A) = e'-""^-'^ lim lE-Je" 

A^^O A^^O L 



A^.^() L 



= IE 



[exp{v(C"-C°)+iACi^} 



□ 



We are now able to end the proof of Theorem 12.21 (item 2). We may apply, without any 
change, the measure tension criterium used in the proof of convergence of (Zf^) in the case 
ao ~ Po = 1 (see the end of Section|4]). This, and Proposition 16 . 61 show that (5^)t>o converges 
in distribution as A^; — > to the Brownian motion with drift {£,t)t>o- 

6.3 Proof of Proposition [2731 

We suppose ao — /3o — 1, ci — co < and rAt — A^ where r > 0. 

We briefly sketch the proof of Proposition 12.31 The approach is similar to the one developed 
in the case 2) of Theorem 12.21 We only prove that converges to the Gaussian distribution 
with 0-mean and variance equals —rcot. Using Theorem 3 in [4], it is equivalent to show 



^lim^E-i [e""^' J = e ! — , V/i G : 
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We have already observed that we may reduce to the case t/At G IN; in this case we have 
= Zt and 

where <^{\,t) is the moment generating function associated with (Xt) (see the beginning of 
subsection 16. If) . Recall that $(A, t) is given by identity (|6.15|) . 
Note that: 

Q = QO + CqA^ = 1 + CqAi, (3 = (io+ CoA^ = 1 + CqAj;. 

Since ct and /3 have to belong to [0, 1], this implies that co < 0. Recall that T>, ■&+ and ?9_ are 
the real numbers which have been defined by (|6.12|) resp. ((6A4)) (with A = e""^-). We have: 

V = 4cgA^ cosh^(/iA,) + 4(1 + 2coA:,), 

i9± = -coA^ cosh(AtA^) ± y^cgAi cosh^(/iA^) + l + 2coA^. 
Using classical analysis we get: 

yO/2 = yiT2^^jArT^^AiT^^(Ai) = 1 + coA, + o(A^), 

^+ = 1 - ^A^ + o(A^), ^_ = -1 - 2coA, + o(A,). 

lim a_|. = 1, lim i?'^^* = lim exp [ ^ A^^l = exp I — car— A, 

lim a_=0, lim \-d-f^^^ lim exp ( 2co A^ 1 = lim exp(^^].=0 (co < 0). 

A^^O A^-.o' ' A^^O I At J A^^O L Ai J 

Relation (|6.15p implies that the variable Z^ is asymptotically normal distributed with 
variance —rcot. 
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